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Abstract 

We investigate the effective potential of the PT symmetric \—g<j)\ field theory, perturbatively 
as well as non-perturbatively. For the perturbative calculations, we first use normal ordering to 

O 

obtain the first order effective potential from which the predicted vacuum condensate vanishes 
exponentially as G — » G + in agreement with previous calculations. For the higher orders, we 
employed the invariance of the bare parameters under the change of the mass scale t to fix the 
transformed form totally equivalent to the original theory. The form so obtained up to G 3 is 
new and shows that all the 1PI amplitudes are perurbative for both G -C 1 and G> 1 regions. 
For the intermediate region, we modified the fractal self-similar resummation method to have 
a unique resummation formula for all G values. This unique formula is necessary because the 
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effective potential is the generating functional for all the 1PI amplitudes which can be obtained 
via d n E/db n and thus we can obtain an analytic calculation for the 1PI amplitudes. Again, the 
resummed from of the effective potential is new and interpolates the effective potential between 



the perturbative regions. Moreover, the resummed effective potential agrees in spirit of previous 
calculation concerning bound states. 
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I. INTRODUCTION 



In recent years, it has been established tha ; 



turn models have real and discrete spectra 



;he PT symmetric and non-Hermitian quan- 
ta 
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This may draw the attention to the reinvestigation of the previously rejected non-Hermitian 
models especially for the quantum field versions. For those theories, there are some clues 
concerning the importance of the study of them. For instance, the very recent result con- 
cerning the existence of a non-Hermitian representation for Hermitian theory which has 
simpler calculation II q . Accordingly, one may aim to describe hadrons in a simple way 
rather than the complicated QCD model by finding a non-Hermitian representation for 
QCD. However, this is still a speculation and in this work we concentrate on the study of 
the effective potential of a simple model but exhibits non-trivial features. A class of simple 
but non-trivial quantum mechanical models are given by 

H = p 2 + x 2 (ixf , e>0. (1) 

All such models have real and positive spectra even in the case of e = 2. In fact, all the 
complex PT symmetric Hamiltonians have real and positive spectra Q]. 

Unlike the quantum mechanical versions of the PT symmetric models, the quantum field 
versions of such theories still are not investigated in a sizable way. Moreover, the quantum 
field Hamiltonians have interesting properties. For instance, a simple model with Lagrangian 
density like 

L = \ ((d0) 2 - m 2 2 ) + f 4 , (2) 

exhibits asymptotic freedom as well as having bound states Q|. Also, for certain range 
of the coupling values it has two-body bound state (meson-like) and for another range it 
has three-body bound state (baryon-like). Thus, it is concluded that this simple model 
has supersymmetric features. This may give us a hope to describe strong interactions with 
Abelian theories without the need for Glouns. Regardless of these legitimate hopes, the PT 
symmetric and non-Hermitian theories, like any of the physically acceptable models, deserve 
the employment of the usual machinery of investigation under different conditions. For 
example, we need to know how they behave at zero and non-zero temperatures, the presence 
of external sources, extreme conditions, etc... . Since the effective potential serves as the 
generating functional for all the one-particle irreducible (1PI) amplitudes, it's investigation 



is the basic stone for all other discussions. Up to the best of our knowledge, the effective 
potential of the non-Hermitian has never been obtained in a form reliable in all regions of 
the coupling space. This situation is partly due to being a new field of study and partly due to 
the non-Borel summability of the theory because of the existence of classical soliton solutions. 
In this paper, we offer a coherent formula for the effective potential of the PT symmetric 
and non Hermitian 4>f +1 theory which is reliable for any coupling value. First, we study the 
effective potential (at zero temperature) of the model in Eq. (J2J) in an effective quasi-particle 
theory which verify perturbation for both g <C 1 and g>l. For the intermediate region in 
the coupling space, the quasi-particle theory ought to be non-perturbative and to resum the 
perturbation series one has to resort to a resummation technique rather than Borel technique 
because the theory is not Borel summable. Pade approximation are suggested for the sake 
of getting reliable results from the input information of perturbation series. However, the 
knowledge of only a few first terms does not permit one to use these techniques jisj]. This 
means that to provide a reasonable accuracy, these techniques need to know tens of first 
terms of the perturbation theory Q|. In fact, going to higher orders in quantum field 
models is not an easy task as time-ordering of many fields results in many different types 
of Feynman diagrams and thus one needs (if it is possible to do the calculations) a long 
time to accomplish the diagrams calculations. To overcome such difficulties, the self-similar 
rod was suggested as a non-perturbative tool for the resummation of divergent series 
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22,123, 



24. 



25L l2fi| . Although this method can give good results even with few 
terms of perturbation series, sometimes the method is not applicable at all. We will argue, 
later in this work, it's applicability and suggestions for modification to render it applicable 
for the effective potential for any coupling value. 

The paper is organized as follows. In Section [HI we obtain the first order effective 
potential by normal ordering of the fields with a field as well as mass shift. Then, the first 
order calculations is supplemented by perturbative corrections up to g 3 . In Section ITTT| we 
briefly review the key points of the fractal self-similar method and introduce more control 
functions to render it applicable for the effective potential for all the values of the coupling 
g. In Section llVl we present and discuss the results and Conclusion follows in SectiorfVl 
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II. THE PERTURB ATIVE EFFECTIVE POTENTIAL 



In low dimensional super-renormalizable theories, it is often enough to work with normal 
ordering to render the quantum field theory finite. This is because there are only few 
diagrams that are divergent and these are regulated by normal ordering. The (^<^ 4 ) 1+1 
theory is such an example that has only one divergent diagram in the self-energy amplitude. 
In that case, one shall start with a Hamiltonian that is normal ordered with respect to the 
vacuum of mass parameter m. 

H = N m {^- ((V0) 2 + vr 2 + m 2 2 ) - f 4 ) . (3) 



We can use the relation [2 



N m exp (i(34>) = exp (~^ 2 a) N M=vl . m exp {if3<j>) , (4) 

to rewrite the Hamiltonian normal ordered with respect to a new mass parameter M = y/t-m. 
In eq.(HJ), expanding both sides and equating the coefficients of the same power in (3 yields 
the result 



N„ 



N. 



N m <p 2 = N 2 M( p + A, 
N m <p 3 = N M <f> 3 + 3AN M c 



(5) 



N 



+ QAN M (p 2 + 3A 2 



with 



A = Int. 

47T 



(6) 



Also, it is easy to obtain the result 



-7T 



N 



M 



1 (V0) 2 + \^ + i- (M 2 - m 2 ) . 



2 v ' 2 J \2 
The mass shift m — > M should be accompanied by the canonical transformation 



(7) 



28] 



(05 7r) — > (-0 + B, IT) 



(8) 
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The field ip has mass M = ■sjt ■ m , B is a constant, the field condensate and II is the 
conjugate momentum (-0). Therefore, the Hamiltonian in Eq.Q can be written in the form; 



H = H + H I + H l + E, (9) 

where 

Ho = N M Q (n 2 + (VVO 2 )) + \n m (m 2 - 3g {B 2 + A)) ^\ 
H I = ^-N M (^ + 4B^ 3 ). 
Hi can be found as 

H x = N M (m 2 - g (B 2 + 3A)) B^, (10) 
and the field independent terms can be regrouped as 

1 / 2 12#A\ 2 4 1 , 2 , 1 2a /-.-.N 

£/ = - \ m B B H M — m ) — m A. 11 

2V 4/ 4 8vr v ; 4 2 v; 

Taking b 2 = AnB 2 and the dimensionless parameters t = G = 2 ^ m2 , the corresponding 
vacuum energy density can be written as 

E(b, t,G) = ^ U 2 - j (b 4 - 6b 2 In t + 3 In 2 t) + t - 1 - In tj (12) 

The renormalization conditions are given by 

—E(b,t,G)=g n , (13) 
where g n is the ip n coupling. For instance, 

pp. rp 

where #2 — ^ 2 — m2 — 3g (I? 2 + A). Note that, the renormalization condition §§ = 
enforces H\ to be zero. 

The quasi-particle Hamiltonian in Eq.flHJ) exhibits some interesting properties. For in- 
stance, the renormalization conditions predicts an imaginary condensate which turns the 
Hamiltonian to be non-Hermitian. To show how this comes out, consider the equations 

BE 

— = (2 + (-G)(b 2 -3\nt))b = 0, (15) 
^ = (_G)(36 2 -31nt)+2 = 2t. (16) 



For b ^ 0, they simplify to 



2- G{\? - 3 hit) = 0, 
-G(36 2 -31nt) +2 = 2t. (17) 



Now, Eq. (|T7jl can be parameterized as 



31nt' 

b 2 = -~ (18) 

This parameterization shows that the parameter t should be less than one and b is pure 
imaginary. This result clearly shows the non-Hermitian property of the theory considering 
the form of Hi in Eq.Q. 

The solutions of Eq. (JT7|) are given by 

b\G) = -3W , (19) 

t(G) = 3GW (^e-&) , (20) 



3G 



where the 
expansion 



ambert's W(x) function is defined by W{x)e w ^ = x. Also, W has the series 
ion 3(| 

H , W = E H^£, (21) 



n>l 



which is convergent if \z\ < \- Thus, for G — >• + , Eq.(fIU)l takes the form 

b\G) ~ -3 (i e w) , (22) 
6 = ±iV^(-^=e^V (23) 



This result predicts pure imaginary condensates which vanishes exponentially as G — > 



17|. Another interesting property that the 



which agrees with the prediction of Refs. Q, 
quasi-particle Hamiltonian posses is that it is totally equivalent to the original theory in 
the sense that setting the parameter t equal to one (M = m), the Hamiltonian H (Eq.®) 
reproduces the original form in Eq.fjSJ). In fact, this is very important because, as we will 



see later in the work, the direct calculation of higher orders, spoils out this equivalence. 
Also, since normal ordering can not account for non-cactus Feynman diagrams, we resort to 
renormalization group invariance to fix the parameters in the theory which automatically 
turn the quasi-particle Hamiltonian equivalent to the original Hamiltonian. Besides, it 
is interesting to note that Eq.([17j) can be obtained by adding a counter term as well as 



22]. 



employing the invariance of the bare couplings on the change of the mass scale t 

In spite of all of the above correct features, the normal ordered effective potential in 
Eq.(fT2j). as we will discuss, is non-perturbative for intermediate values of the coupling G. 
In order to improve the representation of the effective potential near the non-perturbative 
region, we consider the modification of Eq.(|12|) resulting from the higher order perturbative 
corrections to the vacuum energy followed by a modified fractal self-similar method (the 
theory is not Borel summable). 

The normal ordered effective potential of |0 4 theory (Eq.{|12|)) agrees with GEP results 
which in turn accounts not only for the leading order diagrams but also for all the non- 
cactus diagrams 0,0]. Thus, to go to higher orders we include only non-cactus diagrams 

(Fignj). 




FIG. 1: The vacuum diagrams (up to g 3 ) of the effective quasi-particle (— f </> 4 — gBcf) 3 ) theory. 

In the equivalent quasi-field theory, the interaction term is (— 10 4 — g50 3 ). Up to g 3 , 
we have the Feynman diagrams (non-cactus) shown in Fig^ Accordingly, the perturbation 
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corrections to the effective potential are 

8nE(b,t,G) , j2 1 ^ /1, 4 3, 2 3, 2l v 

t-lnt + 6 2 -l-G -6 4 + -ln 2 t- -fr 2 lni (24) 



/;/- " V I ' 4" 2 

+ £? ,_ 3.515^ -0-^ + 9.918* 



t t J \ t 2 t 2 

In fact, this form of effective potential does not predict the values of b and t parameters that 
makes the quasi-particle theory equivalent to the original theory. To kee p th e equivalence, 
we use the fact that the bare parameters are independent of the scale t j36|. Accordingly, 
we obtain the result 

8nE(t,b,G) , l2 „ „ fl lA 3 , 2 3, 2l 

v ' ' - t-\nt + b 2 - 1- G i -b 4 + -\n 2 t- -b 2 \nt 



m 2 V 4 1 2 

+ G 2 (-3.155 3.5156 2 I ^ - 1 } ) ( 20 ) 



- G 3 (^4.057 ^_ - lj + 9.9186 2 ^ - 1 

We employ the conditions in Eq.()13|) to predict the parameters b and t which define the 
vacuum. This leads to the following conditions on the parameters b and t (for b ^ 0) : 

2 - 7. 03G 2 Q - lj - 19.836G 3 - lj -G(b 2 -3\nt) =0 (26) 
-19. 836G 3 f i - lj - 7. 03G 2 Q- - 1 j - 3G (b 2 - hit) + 2 = 2t (27) 

In all of the G 3 order equations (Eq.()25p. Eq.(j2HJ) and Eq.(|27|)) setting t = 1, we have 
b = E = 0, which when substituted back into Eq.Q we get the original Hamiltonian form in 
Eq.©. Up to the best of our knowledge, such form of the effective potential which conserve 
the equivalence between original theory and the transformed one has never been obtained 
before (up to g 3 order). 

As in Eq. ()18|) . Eq.(ffi)|) and Eq. (}2Tj) can be parametrized as 



t = -3G\nt + 7. 03G 2 Q - lj - 19. 836(-G) 3 (J^ - lj - 2, 

b 2 = ~ (28) 

This parametrization is exactly the same as in Eq. ()18|) except that the t parameter takes it's 
new form according to the radiative corrections. One may conjecture, relying on this result, 
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that the prediction of the vacuum condensate as b 2 = ^ is correct up to any order of the 
perturbation series and the higher orders just redefine the parameters due to the radiative 
corrections. 

As we will see in section IIV| the G 3 order equations so obtained has an interes ting 
property. These equations are perturbative not only for G< 1, which is proved in Ref.|l7| 
but also for G 3> 1 (see Fig0 Figfjand FigEJ). To obtain an effective potential which is 
reliable for all the coupling values and noting that the theory is not Borel summable, we 
resort to the fractal self-similar method Q|. In fact, we will modify it via the introduction 
of more control functions to make it applicable for all the regions in the coupling space. 



III. THE MODIFIED FRACTAL SELF-SIMILAR METHOD 

To get finite values of a divergent series, the so-called resummation methods are suggested. 
The most often used among such techniques are the Borel summation |37l LlSl and the 

n n nn 

construction of Pade approximants [39], including the two-point |4£| and multivalued |41l. I42| 
Pade approximants. These techniques have many known limitations. For instance, The Borel 
method can not be applied in case of the existence of classical soliton solutions [31] like the 
model we study in this work. Also, Borel and Pade techniques require to have a number of 
perturbative terms which often are hard to get. Rather than that, the self-similar method 
can give good approximation with few terms at hand. One of the basic ideas in self-similar 
approximation theory is the introduction of control functions which govern the evolution 
of an approximation dynamical system to be close to a fixed point. To introduce control 
functions into a given asymptotic series, one has to employ a transformation that include 
trial parameters. This transformation has to simulate the self-similarity property hidden in 
the given perturbative sequence. For power series, it looks natural to emp loy the power- 
law transformations Since power laws are typical of fractals 0, Q| the power-law 



transformation can also be called the fractal transformation |43j. Accordingly, it may be 
more plausible to introduce control functions via a fractal transform because it satisfies the 
scaling relation 

P{Xx,s) _ xs P(x,s) 



p(Xx) p(Xx) 

which is a typical of fractals, where P(x, s) and p(x) are the polynomials defined below. 
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For a function f(x), the fractal transform is given by 

F(x,s)=x s f(x), (29) 

and it's inverse transform is: 

f(x) = x- s F(x,s). (30) 

Here we follow the work in Ref . to present the key points of the self-similar method. 
Consider the series given by 

k 

Pk = ^2a n x n . (31) 

k=0 

Applying the fractal transform we get 

k 

P k (x, s) = x s P k (x) = a n x s+n . (32) 

n=0 

Define the initial approximation Pq (x, s) = ciqx s = f. Solving for x we get 

x(f,s)=(l)\ (33) 

Then 

k / \ n 

Vk (/, s) = P k (x (f, s) , s) = V a n ( -f- ) S+1 . (34) 

Note that, self similarity means that 

Vk+p (/, s) = y k (y p (f, s),s). (35) 

Considering the cascade y k as a dynamical system with the time as k, then the cascade 
velocity is given by 



k f f \ ^~ ^ f f \ 

ilk (/, s) - y k -i (/, s) = "V a n ( — J - Y] a n I — J 

\ a o/ „_ n \ a o/ 



n=0 v u/ n=0 



a fc ( — ] ■ (36) 



After introducing the control functions, the regime of the self-similar renormalization is 
to consider the the passage from one approximation to another as a motion with respect to 
the approximation number k = 0, 1, 2, .. .. In fact, the trajectory y k (f, s) of this dynamical 
system is bijective, that is, in one-to-one correspondence to the approximation sequence 
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Pk (x, s). This dynamical system with discrete time k is called the approximation cascade. 
The attracting fixed point of the cascade trajectory is, by construction, bijective to the limit 
of the approximation sequence Pk (x, s), that is, it corresponds to the sought function. 

One can deal with continuous time t rather than the discrete time k such that the trajec- 
tory passes through the same points when t = k. In this case, the flow velocity is governed 
by a differential equation rather the difference equation in Eq. (}3l)j) . In other words, the 
evolution equation for the flow reads 




(37) 



Accordingly, the evolution integral is 




(38) 



Thus, the self-similar approximation is given by |4j 



\ sa / 

where t* k = 1 when no restrictions are imposed on the series 
The applicability of the method is governed by the stabilizers 




(39) 




(40) 



or their images 




(41) 



The stability condition is given by 




(42) 



For the series given above we have 



k 




(43) 



For k = 3, the stabilizers are given by 



m k (x, s) 



xa\ + 2x 2 d2 + 3x 3 a3 ao + xa\ + x 2 ci2 + x 3 a^ 



(44) 
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The most stable aproximant is obtained if m k (x, s) = 0, or 

Mi + 2x 2 a 2 + 3x 3 a 3 



(45) 



do + xa\ + X 2 Cl2 + X 3 d3 

If m k (x, s) does not have a positive root, then m k (x, s) is a monotonic decreasing function 
of s. Therefor, the minimum is given by 



\m k (x,s)\ 



cio + xai + x 2 a 2 + x 3 a3 



«o 



(46) 



If all \m k (x, s)\ are less than one, the bootstrap formula given by 



a>i a 2 f a 3 f a 4 f a s f a k \\\\\\ 
p k = a exp — xexp — xexp — xexp — xexp — x exp x (47) 

\ a o \ a i V«2 V a 3 V a 4 \a k -i J J J J J J 

represents the resummed series j^J. However, it is not guaranteed to have all \m k (x,s)\ 
less than one for s — > oo and for every point of the argument x. In fact, this is the 
situation in case of applying the fractal self-similar method to the effective potential in 
Eq.(|25J) where there are some coupling values for which the bootstrap formula does not 
exist. So, instead of applying the the fractal self-similar method directly to the series in 
Eq.(|25J) we introduce another control function via a bijective transformation which has the 
property of transforming the original series to another one for which the bootstrap formula 
is applicable. Then, we apply the fractal self-similar method to the transformed series and 
at the end we apply the inverse of the transformation to get the resummation formula of 
the original series. 

To test the modification we introduced, consider the Lambert W function defined by 

W(x)exp(W{x))=x. (48) 
The series expansion of W (1 + x) is 

W (1 + x) » W (1) + - ~~~-T~r 3? +(-l(W (l)) 2 2 + W 0) \ x * (49) 



1 + W(1) V 2 V v 77 (1 + ^(1)) 
+ ( V (I)) 3 9 + 8 ^ (1) + 2( T (1)) 1 - 3 + O W . (50) 

At x = 3, W (1 + x) = 1. 202 2 and the perturbative result (up x 3 ) is 1. 918 9. The error 
percent is | L2022 -^ 2 9189 1 % = 59.6 16 %. 
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Let us apply the transformation T (W (1 + x)) = W (1 + x) + c, where c is used as a 
control function. Apply the fractal self-similar method to T (W (1 + x)) and find c which 
make all the \rrik (x, s)\ less than one and then apply T -1 to the obtained bootstrap 
formula we get the result W (1 + x) ~ 1. 1798 with the error percent | L2 °i 202 2* 798 I 100% = 
2. 069 7%. This result indicates the success of the modification we introduced to the fractal 
self-self-similar method which can be summarized as: instead of using one control function 
s we introduce two control functions s and c. With s — > 00, we adopt c to obtain stable 
approximant. In fact, this trick is necessary because it results in a unique formula for the 
approximant for all the values of the coupling G. Since it is well known that the effective 
potential is the generating functional for all 1PI amplitudes, unique formula for the effective 
potential makes it easy to obtain the different amplitudes via analytic differentiation. 

We applied the modified fractal-self similar method to the perturbation series of the 
effective potential in Eq.()25|). As we will see in the following section, the resummed series 
fits the perturbative data for regions where the perturbation series is reliable (G < 1 and 
G> 1). 

IV. NUMERICAL CALCULATIONS AND DISCUSSIONS 

In this section, we present the numerical calculations concerning the perturbative as well 
as non-perturbative calculations for the vacuum properties of the (^0 4 ) 1+1 non-Hermitian 
field theory. The non-Hermiticity of the theory is clear in it's quantum field version without 
the employment of boundary conditions in the complex x space. This can be extracted from 
Eq. f)28j) , where we realize that the condensation is always pure imaginary and thus the from 
in Eq.© shows that the theory is non-Hermitian as well as PT symmetric. 

In Fig|2J the effective potential is plotted as a function of the coupling G for different 
values of order of perturbation k as well as the modified fractal-self similar resummation 
formula. Careful analysis of the plot shows that the theory is really perturbative for both 
G < 1 and G ^> 1. Although we can extract visually from the graph that the non- 
perturbative calculations coincides for k = 1,2 and 3 for G <C 1, it is not clear that the 
theory is perturbative for G ^> 1 region. To clarify the prturbative character of the theory 
for G ^> 1, we rearrange Eq. pp} as: 
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m 

where 



^^ + * , (51) 



-3.5156 2 



G 9.9186 2 



1 



(52) 



E = t - In t - 1 - G (j In 2 ^ + G 2 ^-3.155 Q - 1^ - 4.057G 3 (J^ - l^j . (53) 

In fact, the field dependent term Ej, behaves well for both G <C 1 and G> 1. For G 1 
it goes to zero as G — > while £7q goes to infinity. The appearance of infinite effective 
potential as G — > can not be considered as Infra- Red divergence because it appears only in 
the vacuum energy which can be rescaled as we can measure only differences in energies. All 
the 1PI amplitudes are finite as G —>■ because we can get them by successive differentiation 

n 

of the effective potential with respect to the condensate b |29f while E has no contribution. 
Thus, we can safely consider the effective potential as Eb only. Also, when G ^> 1 the 
theory is perturbative (see FigOJand FigEJ). 

The fractal-self similar resummation for Eb is plotted in FigJ3] where we can realize that 
it goes to zero as G — > and decreases as G be very large. In fact, this is an indication that 



as we increase the coup 
proved before in Ref. 



ing G the number of bound states decreases as G increases which is 
17| . To assure this point, we plot the expectation value of the potential 
term in FigJBl which shows that the depth of the potential decreases for large G values and 
thus it is expected to have no bound states for G ^> 1. 

V. CONCLUSION 

We investigated the effective potential of the (— 10 4 ) 1+1 in view of the effective potential 
representation. The normal ordered effective potential predicts a pure imaginary condensate 
but tends to zero (exponential decrease) as g — > + . The imaginary value of the condensate 
turns the theory non-Hermitian without a priori postulated boundary conditions which is 
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necessary in the quantum version of the theory. The calculation of the effective potential is 
extended perturbatively up to g 3 order which in turn predicts the same shape for the vacuum 
condensate but with higher values. Up to g 3 , we found that the theory is perturbative for 
both j < 1 and g ^> 1 regions. Since the theory is ought to be non-perturbatieve for 
intermediate regions of the coupling space, we supplemented the perturbation result by the 
fractal-self similar method in a way that results in a coherent formula to make the calculation 
of all the 1PI amplitudes accessible analytically. 

We believe that these new results concerning the effective potential of the PT symmetric 
4>f + i theory is a step forward toward the understanding of these novel kind of field theories. 

A very interesting note should be mentioned, the simple model we used shows many 
similar features with QCD. For instance, In Ref.^j, it was shown that the theory has 
two body bound state similar to qq and three body bound state like baryons. Also, our 
calculations shows that the theory is perturbative for large scales (Fig0J) which may be a 
clue that the theory may have an Ultra- Violet fixed point like QCD. Finally, the theory 
may have a symmetry restoration for large G values. All of these expected features will be 
our task of the next work. 
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FIG. 2: The effective potential 8tt^ up to G 3 as a function of the coupling 




FIG. 3: The vacuum condensate squared versus the coupling G for first, second and third order in 
the perturbation series. 
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FIG. 4: The reciprocal of the 1PI two-point function versus the coupling G for first, second and 
third order in the perturbation series. 



E -0.6- 



LU 

CO 




FIG. 5: The resummed vacuum energy Ej, as a function of the coupling G. 
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FIG. 6: The expectation value of the potential term in the Hamiltonian as a function of the coupling 
G. 
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